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^ ' Abstract 

^ ■ Measurements of polarization fractions in B ^ VV transitions, with V a light 

^ . vector meson, show that the longitudinal amplitude dominates in 5° -^ P^P^ , 

H ! B^ — > p'^p^, and B^ -^ p^K*^ decays and not in the penguin induced decays 

■ " " ' 5° — >• (pK*^, B^ -^ (pK*^ . We study the effect of rescattering mediated by 

charmed resonances, finding that in i? — >■ (pK* it can be responsible of the 

suppression of the longitudinal amplitude. For the decay B —>■ pK* we find 

that the longitudinal fraction cannot be too large without invoking new effects. 



1 Introduction 



An important result obtained by Belle and BaBar Collaborations is the measurement of 
the decay widths and of the polarization fractions of several B decays to two light vector 
mesons [1, 2, 3, 4]. The branching fractions measured by the two Collaborations are 
collected in Table 1 together with the averages. Together with these data one should 
collect the upper bound B{B^ -^ P^P^) ^ 2.1 x 10~^ from BaBar [2]. Through the 
analysis of angular distributions, the polarization fractions of the final states have been 
measured as reported in Table 2. In the decay modes B^ -^ P^P~ and B~^ -^ P^P~^, p^K*^ 
the final states are essentially in longitudinal configuration, with a larger uncertainty for 
B^ -^ p^K*^\ on the contrary, in both the observed B -^ (pK* transitions the longitudinal 
amplitude does not dominate, providing nearly 50% of the rate. 



Mode 


Belle [1 


BaBar [2 


Average 


B+ -^ (t)K*+ 


{Q.7tti^%) X 10-« 

iw-o^mi) X 10-6 


(12.71^:^ ±1.1) X 10-6 
(11.2 ±1.3 ±0.8) X 10-6 


(9.5 ±1.7) X 10-6 
(10.7 ±1.2) X 10-6 


Mode 


Belle [3 


BaBar [2, 4 


Average 


B+ ^ p^K*+ 
B+ -^ p^p+ 
B^ -. p+p- 


(31.7±7.ll|?) X 10-'^ 


(10.61^:^ ±2.4) X 10-6 

(22.5+^i±5.8) X 10-6 

(25™) X 10-6 


(26.2 ±6.2) X 10-6 



Table 1: Branching fractions of -B ^ VV decay modes. 



There are reasons to expect that the light VV final state should be mainly longi- 
tudinally polarized, see, e.g., the discussion in [6]. In the following we summarize the 
arguments, which essentially rely on factorization and on the infinite heavy quark mass 
limit. Invoking such arguments, the deviation observed in B -^ (pK* could be interpreted 
as a signal of new physics [7]. A more orthodox interpretation [6], in the framework of 
QCD improved factorization [8], relies on the observation that (logarithmically divergent) 
annihilation diagrams can modify the polarization amplitudes in B ^ 4>K*, producing 
fractions in agreement with observation. 

In this note we wish to address another effect that potentially changes the result in 
the penguin induced B -^ (pK* decay without affecting the observed B ^ pp transition: 
rescattering of intermediate charm states. Such effects, studied long ago in B ^ Kir 
transitions [9] and investigated recently in other B -^ PP and VP decays [10] as well 



as in factorization forbidden B transitions to charmonium final states [11], can invalidate 
the arguments on the basis of which the dominance of the longitudinal configuration is 
argued. 

We discuss factorization and its consequences in Section 2 and the analysis of rescat- 
tering effects for B^ -^ (pK*^ in Section 3. At the end we discuss a few consequences. 



Mode 



Pol. fraction 



Belle [1] 



BaBar [2] 



Average 



fi^ 



(l)K*~ 



5° -^ (pK*^ 



Vl/V 0.46 ±0.12 ±0.03 

Vl/T 0.43 ±0.09 ±0.04 0.65 ± 0.07 ± 0.02 0.58 ± 0.06 

(0.52 ±0.07 ±0.02) 

r^/r 0.41 ±0.10 ±0.02 (0.27 ±o.o7± 0.02) 



Mode 



Pol. fraction 



p+p 



Belle [3] 



BaBar [2, 4] 



Average 



0.96"'^± 0.04 
0.97to:o7 ± 0.04 
0.98to:ol ± 0.03 



B+- 
B+ 
B' 



Tl/F 



0.95 ±0.11 ±0.02 



0.96 ±0.07 



Table 2: Polarization fractions in i? — > VV transitions. The BaBar results reported in 
brakets are preliminary data quoted in ref. [5]. 



2 Polarization in factorization-based approaches 

The decay B^ -^ (j)K*^ is described by the amplitude 

A{B\p) ^ <P{q, e)K*\p', r^)) =A,e*-r^*+ A (e* ■ p) {v* ■ q) + Mi e^^^' Cr^^p^p', (1) 

with e(g, A) and ri{p' , A) the and K* polarization vectors, respectively, with A = 0, ±1 
the three helicities. Since the decaying B meson is spinless, the final vector mesons share 
the same helicity. ^o and A2 are associated to the S- and D-wave decay, respectively, and 
Ai to the P-wave transition. 

The three helicity amplitudes Al and A± can be written in terms of .4.0,1,2: 

1 



Al 



M^Mk 
-AqTMbW\Ai ■ 



[{p-p'-Ml,)A, + MJ,\f^\'A2 



(2) 



in the transversity basis, the transverse amphtudes 

A = ^^^ = -V2A, (3) 

A^ = ^^^- =-V2MBW\Ar 

can also be defined, with \p\ = A2 (M^, M|^*, M|)/2Mb (A the triangular function) the 
common and K* three-momentum in the rest frame of the decaying B-meson. In terms 
of such amplitudes the expression of the decay rate is simply: 



SnMl 



(1^,1^ + 1^1,1^ + |^J^) , (4) 



r 

rii 


\Al 


1^11 1' 


' + \A^' 


r 


Al' 


A^\ 


+ \A^ ' 
2 



while the three polarization fractions are given by 

h - 

^^ ~~ r \Al\^ + \Aii\^ + \a±\^ ■ 

In order to compute the amplitude eq.(l), we consider the effective weak Hamiltonian 
inducing the b — > sss transitions, which can be written as 

/-f 10 

Hw = -^{-V*Vt,){Y.c,0, + c^Or, + Cs,Osg) (6) 

V^ i=3 

with the operators 

q' 

q' 

O5 = {baSa)v-A^{qi3q'i3)v+A 
<?' 

Ce = {bpSa)v-Aj2i^a(l'p)v+A (7) 

q' 

3 

^7 = 2^KSa)v~Aj2^q'(^^'i'(s)v+A 

q' 
3 

^8 = -ib(3Sa)v-Aj2^q'(9W^)v+A 

q' 



3 

Oq = ^{baSa)v-Aj2^Q'(^f3Q'f3)v-A 

3 

q' 

{a,P are colour indices and {qq)vTA = Ql'^i^ T 75)q')- C^3-6 are gluon penguin oper- 
ators, O7-10 electroweak penguin operators, Oj^ = -^mjja^^ {1 + 75)sF^,y and O^g = 
-^rribba^^ (1 + 75)T"sG'^^, with F^^ and G^^, the electromagnetic and the gluon field 
strength, respectively; Cij^^sgif^) are the Wilson coefficients. 

The amplitude A{B^ — > (pK*^) obtained from (6) admits a factorized form 

AfactiB' - 0ir*°) = ^i-V;,Vt.)a^{K*'ip\r^)\ibs)v^A\B\p)){<Piq,e)\irss)v\0) (8) 

with aw = 03 + 04 + 05 — -{a-i + ag + aio), Oj = q + ^— for i = 3, 5, 7, 9 and Cj = Cj + ^— 
for i = 4, 10 (A'^c is the number of colours). This formula presents the drawbacks of naive 
factorization, namely there is not a compensation of the scale dependence between Wilson 
coefficients and operator matrix elements. However, it allows us to immediately write 
down the polarization fractions, once the B -^ K* matrix element has been expressed in 
terms of form factors ^ , and the meson leptonic constant has been introduced: 

{<Piq,e)\s^^s\0) = UM^e*>' (9) 



{K*{p',r^)\h,il-^,)s\Bip)) = -^6^.p.r^>V %, , ,, - [(M^ + M^O^i^ 



2V 

MbTMk* 



^' -{V*-P){p + P'),-2MK-S-^^^{V*-P)q,] 



Mb + Mk*' -^ ^- ' - ^^ - g2 

(10) 

with the form factors V^ Ai, A2, A3 and Aq functions of q^ . From (8-10) it is easy to write 
down the polarization amplitudes and check that, for large values of Mb-, 

Al oc Ml[{A,{Ml)-A,{Ml)) + ^{A,{Ml) + A,{Ml))] 

A\\ oc MbA^{MI) (11) 

Ai. oc MbV{MI) , 

expressions which determine the behaviour of the three amplitudes once the parametric 
dependence on the heavy quark mass of the form factors close to the maximum recoil 



^For the B -^ K* and B -^ D* matrix elements eqs.(lO) and (18) we use the same phase convention. 

5 



point has been established. In the hmit Mb -^ oo and for g^ = such a dependence has 
been investigated [12] with the resuh that the three form factors V, Ai and A2 should be 
equal: A2/A1 = V/Ai = 1. One therefore expects: 

regardless, in this scheme, of the Wilson and CKM coefficients. Assuming generalized 
factorization, with the substitution of the Wilson coefficients Oj with effective parameters 
a^ , it is eventually possible to reconcile the branching ratio with the experimental mea- 
surement, but not to modify the polarization fractions, since the dependence on the a, 
cancels out in the ratios. Therefore, in order to explain the small ratio F^/r within the 
Standard Model one has to look either at the finite mass corrections, or at effects beyond 
factorization. 

For finite heavy quark mass, one can compare the experimental result for the polariza- 
tion fractions in B^ -^ (f)K*^ decays (Table 2) with the predictions of various form factor 
models [13, 14, 15, 16]. As shown in fig. 1, in many models the ratios A2/A1 and V/Ai 
deviate from the asymptotic prediction, suggesting that the regime of finite Mb does not 
fully coincide with the asymptotic regime. In one case there is a marginal agreement 
between the form factor model and data. However, the indication of effects beyond naive 
and generalized factorization is clear. 

3 Rescattering effects 

If one considers the possibility of rescattering effects, there are other terms in the effective 
weak hamiltonian that can induce the transition B^ -^ (j)K*^. Processes that should be 
the most relevant ones are b -^ ccs -^ sss. Such processes can give sizeable contribution 
to the penguin amplitudes obtained from (6) since they involve Wilson coefficients of 0{1) 
(that multiply current-current quark operators), while the Wilson coefficients in penguin 
h -^ sss operators are smaller (0(10^^)). On the other hand, there is not a CKM 
suppression in such processes, since |V^(feVis| and iV^tVcsl are nearly equal. An example of 
processes of this type is depicted in fig. 2, where a sample of intermediate charm mesons 
is shown. As far as the polarization of the final state is concerned, one has to notice that 
different intermediate states in fig. 2 contribute to different polarization amplitudes, so 
that the longitudinal as well as the transverse amplitudes can be modified. For example, 
considering only intermediate pseudoscalar and vector charmed mesons coming from the B 
meson vertex, there are eight diagrams of the kind depicted in fig. (2). Intermediate states 

6 




0.6 0.8 1 1.2 1.4 

A2/A1 



Figure 1: Ratios oi B ^ K* form factors: \/(M|)/Ai(M|) versus A2(M|)/Ai(M|). 
The continuous lines correspond to the (one, two and three-o") regions of the Belle data 
in Table 2 (a) and of the average of Belle and BaBar data (b) for F^/r and F^/F in 
B^ — > (j)K*^. The points correspond to different form factor models: QCDSR [13] (dot), 
LCSR [14] (triangle), MS [15] (square), BSW [16] (diamond). 



comprising one vector and one pseudoscalar meson (four diagrams) only contribute to the 
P— wave transition and therefore to the amplitude A±. On the other hand, intermediate 
states comprising two pseudoscalar mesons (two diagrams) only contribute to Al and A\\, 
while intermediate states with two vector mesons (2 diagrams) contribute to the three 
polarization amplitudes Al, A± and A\\. 

In order to estimate the contribution of diagrams of the type in fig.2 we can use a for- 
malism that accounts for the heavy quark spin-flavour symmetries in hadrons containing 
a single heavy quark [17] and for the so called hidden gauge symmetry to describe their 
interaction with light vector mesons [18]. As well known, in the heavy quark limit, due to 
the decoupling of the heavy quark spin sq from the light degrees of freedom total angular 
momentum Se, it is possible to classify hadrons with a single heavy quark Q in terms of 



B 




Di*) 



4> 



K* 



Figure 2: Rescattering diagrams contributing to B ^ (pK*. The box represents a weak 
vertex, the dots strong couphngs. 



se. Mesons can be collected in doublets the members of which only differ for the relative 
orientation of Sq and Si [17]. The doublets with J^ = (0^, 1^) corresponding to sf = ^ 



can be described by the effective fields 



H„ 



;i+/) 



[PLl' - Pain] 



(13) 



where v is the meson four-velocity and a is a light quark flavour index. The field Ha is 
defined as Ha = '~f^Hl'~f^; all the heavy field operators contain a factor \/Mh and have 
dimension 3/2. 

It is possible to formulate an effective Lagrange density for the low energy interactions 
of heavy mesons with light vector mesons [18]. The interaction term of such a Lagrangian 
reads as 



Chhv = -iP Tr{H,{v^p^),aHa} + I A Tr{H,{a^''' F^,\aHa} . 



(14) 



■gv 



Light vector mesons are included in (14) through the fields p = i — p representing the 
low-lying vector octet: 



(15) 



with F^j, = d^pi, — di,p^ + [p^, py\. Invoking the mixing cjg — ^o one gets the interaction 
term involving 0. The interactions of heavy mesons with the light vector mesons are thus 
governed, in the heavy quark limit, by two couplings (3 and A. From light cone QCD sum 

8 



/ 

1 


[^ rfi -1- 




P+ 




K*+ \ 




P~ 


- 


-^'^p' + \ 




K*^ 


v 


K*- 




K*' 




-Vf^8 ) 



rules [19] as well as from vector mesons dominance arguments [10] one estimates (3 ^ 0.9 
and A ~ 0.56 GeV~^, while gy is fixed to gy = 5.6 by the KSRF relation [20]. 

Using (14) it is easy to work out the matrix elements Dy'D^*'K* appearing in one of 
the vertices in fig. 2: 



{D-{pr>-p')K*\p\r^)\D-{pj, = Mr)VD)) = (3^ MdMdAvd ■ V*) 

{D;ipn-p')K*%p',r])\D*~{pry,Vi)) = tX^M^.M^^ e'^^'^^ VDaVtp'.Vip 
{D:-ipD-p',e^)K*%p',r])\D*~{pD,Vi)) = -P^Md^Md^ (vd ■ vl i^l ■ Vi) 
+yMn*Mn*[ivi-vliel-p') - (vi ■ P) (4 ■ v*)] 

(16) 

where /3 = ^^y and A = — ■^. Matrix elements involving in the other vertex in fig. 2 
are obtained analogously. 

As for the weak amplitude B^ -^ D^A'^ D^*^^ , since there is empirical evidence that 
factorization reproduces the main experimental findings [21], we write it as 

{Di*^+D^*^-\Hw\B-) = ^V,tV:MD^*^'\{V - Ar\B'^){Di*^+\{V - A)^\0) (17) 

with fli ~ 1. In the heavy quark limit the matrix elements in (17) involve the Isgur-Wise 
function [17]: 



{D-{v')\V^'\B\v)) = ^ MbMd i{v ■ v'){v + v'Y 
{D*-{v\t)\V^\B\v)) = -i^MM^,avv')e*pe^^^^vJ^ (18) 

{D*'{v',e)\A>^\B\v)) = ^JMbMd^ i{v ■ v') £^[(1 + v ■ t;')/^ - v^v'^] , 

V and v' being B and D^*' four-velocities, e the D* polarization vector and ^{v ■ v') the 
Isgur-Wise form factor. As for the D^*^ current- vacuum matrix elements defined as 

{Q\qal^^,c\Da{v)) = IdMd^v^ 

mal'c\D:{v,t)) = fD'Mo^y , (19) 

they can be parameterized in the heavy quark limit in terms of a single quantity Jd^ = fc*. 
Now, the estimate of the absorptive part of the rescattering diagrams in fig. 2 

oZtiM^ J-1 

(20) 
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can be carried out. The integration variable z = cos 6' is related to the angle between the 
three-momenta of and of the emitted D'f^ from B vertex in fig. 2. We use | V"cfe| = 0.042, 
\Vcs\ = 0.974 (the central values reported by the Particle Data Group [22]), fo* = fos = 
240MeV[23]ande(t/) = (T^)'. 

The couplings in (16) do not account for the off-shellness of the exchanged D^*^ mesons 
in fig.2. One can introduce form factors: 

9^it)=g^oFit), (21) 

to account for the t-dependence of the couplings (the vertices in rescattering diagrams 
cannot be considered point-like since they do not involve elementary particles), gio being 
the on-shell couplings. However, the form factors are unknown. We use 

m = -j^ (22) 

to satisfy QCD counting rules. We could vary the value of A, considering the uncertainty 
from the form factor F{t) in the final numerical result. Instead, since the relative sign 
of rescattering and factorized amplitude is also unknown, as well as the role of diagrams 
involving excitations and the continuum, we fix A = 2.3 GeV and analyze the sum 

varying the parameter r and approximating the long distance amplitude with eq.(20). 

We compute the short-distance factorized amplitude using the B -^ K* form factors 
appearing in two extreme cases in fig.l, the model [13] and the model [14], with Wilson 
coefficients as = 48 x 10-^ a^ = (-439 - 77i) x 10-^ og = -45 x 10~^ 07 = (-0.5 - 
1.3i) X 10"^, ag = (-94 - 1.3i)48 x 10"^ and a^o = (-14 - OAi) x 10"^ as computed in 
[24] for Nc = 3. 

The result is depicted in fig. 3. For the model [13], a contribution of the rescattering 
amplitude is in order to obtain the measured B -^ (f)K* branching fraction. Of the two 
possible values of the parameter r which reproduce the experimental rate, r ~ 0.08 and 
r ~ —0.3, the former allows us to simultaneously obtain a small longitudinal polarization 
fraction: Tl/T ~ 0.55, compatible with the measurements. The tranverse polarization 
fractions turn out Tji/r ~ 0.30 and r_|_/r ~ 0.15. They are both consistent with mea- 
surement, but with the hierarchy T\\/T > r_L/r. 

If we use the form factors in [14] , for r = the predicted rate exceeds the experimen- 
tal datum, so that the rescattering contribution should be weighted by a negative r to 
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Figure 3: Dependence of the branching ratio and polarization fractions of B^ -^ (f)K*^ 
on the long distance contribution. B -^ K* form factors computed in [13] (left) and in 
[14] (right) are used in the short-distance amplitude, r = corresponds to the absence 
of rescattering. The three curves in (b) correspond to F^/F (continuous curve), F_|_/F 
(dashed) and Fy/F (dot-dashed). The horizontal lines represent the experimental result 
in Table 1 for the branching ratio (a) and for F^/F (b). 



reconcile the branching fraction; as depicted in fig.3, in such a region (r ^ —0.05) the 
longitudinal fraction increases. However, this conclusion crucially depends on the value 
of the Wilson coefficients a^ — aio used as an input in the evaluation of the short-distance 
amplitude. As shown in [24] , for example, 04 varies from —402 — 72i to —511 — 87i 
changing Nc from 2 to 00. For a smaller value of the sum of Wilson coefficients, both the 
sets of form factors would require a similar long-distance contribution, with the effect of 
reducing the longitudinal fraction. 

A feature of both the sets of data is that, in the region of r where the experimental 
rate is reproduced, Fy is larger or similar to F_l. The ratio — ^ is sensitive to operators 
of different chirality which would appear in the effective Hamiltonian in extensions of the 
Standard Model [6]. 
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4 Discussion 

The conclusion of this analysis is that FSI effects can modify the helicity amplitudes in 
penguin dominated processes. The numerical result depends on the interplay between 
Wilson coefficients, form factors and rescattering amplitude, and we have shown that the 
experimental observation can be reproduced. At the same time, the rescattering effects 
we have considered are too small to affect the observed B ^ pp decays. As a matter 
of fact, while the CKM factors in the tree diagram in B^ -^ P^P^ transition {V*ijVud) 
have similar size to the CKM factor in the FSI diagram in fig.2 {V*ijVcd), the Wilson 
coefficient in current-current transition is 0{1). We can expect to observe FSI effects in 
colour-suppressed and other penguin induced B -^ VV decays, such as B^ -^ p^K*^, 
B^ -^ uoK*^, and 5° -^ p^p^, 5° ^ p^uo, 5" -^ p'K*'^, B' -^ K*- K*^ . 

Let us consider B^ -^ p^K*^ . On the basis of general arguments, we cannot assess the 
role of FSI without an explicit calculation, due to the CKM suppression of the factorized 
amplitude. The determination of the rescattering amplitude, similar to that in fig.2, can 
be done following the same method discussed above, obtaining Tl/T ~ 0.7, i.e. smaller 
(even though compatible within 2-0" ) than the measurement in Table 2. 

Therefore, in our approach we can accomodate a small T^ for B — * (pK* at the prize 
of having a smaller value of F^ for B -^ pK*, which is not currently excluded due to 
the uncertainty in the data for this mode. It is interesting to notice that an analogous 
prediction is done in QCD improved factorization [6], where one gets — {B -^ pK*) < 

— {B -^ <pK*). More precise measurements are in order to suggest a solution to this 
polarization riddle. If further measurements of polarization fractions will confirm the 
present situation of a small longitudinal fraction in (pK* and a large longitudinal fraction 
in pK* , in that case we cannot identify uniquely the rescattering mechanism for explaining 
the data, envisaging the exciting necessity of new effects. 
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